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Abstract 

Using the unabridged Maxwell equations (including £, IE, £ and H) new ef- 
fects in collisionless plasmas are uncovered. In a steady state, it is found 
that spatially varying energy density of the electric field (Ej_) orthogonal 
to £ produces electric current leading, under certain conditions, to the re- 
lationship 

Pi * W Ti- =■ 

where e is the dielectric constant of the plasma for fields orthogonal to £. 
In steady state quasi-two-dimensional flows in plasmas, a general relation- 
ship between the components of electric field parallel and perpendicular to 
£ is found. These effects are significant in geophysical and astrophysical 
plasmas. The general conditions for a steady state in a collisionless plasma 
are deduced. With time variations in a plasma, slow compared to ion-gyro- 
period, there is a general current, ( j_*) , which includes the well-known 
polarisation current, given by 

i* = ^ M) + CP X £)] X £ 

where M and £ are the magnetisation and polarisation vectors respectively. 

*0n leave from La Trobe University, Bundoora, Victoria, Australia 3083. 
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Introduction 

In this paper, macroscopic processes in a collisionless plasma are con- 
sidered from the standpoint of Maxwell's unabridged equations including all 
four vectors £ and H. This is in contrast to the widespread approach 

which employs _E and but not £ and H. New results are obtained which are 
not evident in the conventional approach. The paper starts with Maxwell's 
unabridged equations and states, as is well-known, how these are often ab- 
breviated for use in methods which frequently employ velocity distribution 
functions of particles. After uncovering new effects in collisionless plasma 
using general arguments from Maxwell's equations, the approach of this paper 
is placed in context ir plasma physics, and is shown how it extends our under- 
standing of plasmas. 

Maxwell's equations as commonly written are, in c.g.s. units. 


V X E = 


1 

c ¥t ’ 


3D 4ttJ 


V X H = 


c 3t 


(13 

( 2 ) 


where £= electric field, i.e., electric force per unit charge, 

H = magnetic field, 

£ = magnetic induction = H + 47rM, (3) 

M = magnetic dipole moment per unit volume (magnetisation vector) , 

£ = the electric displacement = E + 4 itP, (4) 

£ = electric dipole moment per unit volume (polarisation vector), 

£ = total current excluding displacement current (see Maxwell, 1891). 
Following Maxwell, 


V • D = 4iTp. 


(5) 
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From equation 1 ^ (V * B) = 0, from which it is usual 1> assumed that 

(Stratton, 1941) 

V • B = 0. (6) 

g 

From equation 2 it follows that 47T div J + ^ div !D = 0, and 

7 • J + If- = 0. (7) 


Either one accepts 5 to infer 7 or uses the conservation of charge (equation 
7) to infer 5 (Stratton, 1941), 

Maywell's equations may be presented in other forms (Stratton, 1941) e.g.. 


V X E 

V X B 


1 

"c ’St"' 


, 3E . 3P , 

^3t ^ C^lt 


( 8 ) 

(9) 


V • B = 0 


( 10 ) 


V • E = 4ttp - 4tt V • P 


( 11 ) 


In many plasma physics books and papers (e.g. Spitzer, 1962) instead of 
equation 9, there is written 

3E 


n D 1 477 . 

Vx B = 


C12J 


where, therefore, 


3P 

j = %-jr c V X M + J 

— dt. — — — 


(13) 


and instead of equation 11 is written 


V » E_ = 477Pg, 


(14a) 


3 


wherrs, therefore, 


oRtamf.5- 

OF POOR 


Pe = p “ Z * £ 

for a plasma consisting of electrons (-) and one kind of ion (+) 


(14b) 


j = e (nV - N"V"), 


(15a) 


Pg = CN"" - N“) e, 

V = ^ y"v f (r, V, t) d^v, (15b) 

and 

t) d^v. (16) 

It is argued in texts on plasma physics that the velocity distribution f takes 

into account all the motions of all the particles and hence equation 14a and 

+ - 

integral 14b gives the quantity j_, then using integral 15 to find N and N 
provides the quantity p^. Sometimes £ is taken out of ^ (Clemmow and Dougherty, 
1969) . As expected 


H-f 


f(r, V, 


div j 


3 

3t 


(V • P) + div J 


(17) 



3t 


where p^ = "equivalent” charge density due to polarisation in agreement with 
equation 11. To solve problems of phenomena purely internal to a plasma, 
the equations 12, 13, 14 (a, b) and 15 are commonly used, putting J = 0 and 

p = 0. 
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This paper treats a plasma in quasi -static conditions as a fluid with 
specifiable dielectric constant and magnetic permeability, using the 
unabridged Maxwell equations 1 to 7. The plasma is allowed to be driven from 
the outside by ai)plied currents (J) and charged distributions (p) . This is 
in contrast to the present approach to plasma physics referred to above, 
which, in effect considers the charges and currents to be in a vacuum so 
that £ = ^ and ^ = H (Clemmow and Dougherty, 1969) , and the equivalent infor- 
mation on e and jj is contained in velocity distribution functions and 
appropriate equations of motion. 

Broad consistency between the two approaches is found; however this 
paper brings out new effects driven in a plasma by externally applied charge 
distributions which, in principle, could be quite arbitrary. Suspecting that 
many phenomena in the ionosphere, magnetosphere and interplanetary space are 
driven by outside sources we approach these phenomena theoretically in this 
way in the following discussion. 


Electromagn et ic Stress on a Dielectric Medium 

All discussion from here on is restricted to quasi-static conditions. 

For slow variations of electric field orthogonal to it is known that a 
plasma behaves like a dielectric medium. So that (Stix, 1962} for ^ = (0, 0, 


B ) 
z-' 

D = (eE , eE EJ 
— ^ x’ y, 

CIS) 

where 

47Tp^C^ 

e = 1 + 

C19) 

and 

= i Vs- 

(20) 


n^, m^ = number density, molecular mass respectively of species s of charged 
particle in the plasma. 


5 


ORIGINAL PAGE 13 
OF POOR QUALITY 


Before discus 'iig the macroscopic forces on a plasma let us refer to tlve 
classical discussion of the forces on a volume (V) of fluid with isotropic 
dielectric constant (Stratton, 1941). lie shows that, given J3 = yl[. 


A 

dt 


^^iiech 




f tPB . c”‘ J X B - ^ E"Ve - — ■ I1"7 m 
V 


( 21 ) 


wliere G , is the mechanical momentum of the volume of fluid and G is the 
-mech "<• 

electromagnetic momentum in the volume (V) , The components of the electric 
force ~ ■g™ E^Ve are given, in the case when V x 13 = 0, by 


1 


Z 3 


47TpEj 2 ^ g-- ^ 


1 3 


2 k •' J 


( 22 ) 


In the steady state we may assume from equation 21 tliat 


(4tt)“^ (div 0)b + c"JxB-~ E^Ve 


-1 


1 t*2\ 


= 0. 


(23) 


In addition one can say that the component of ^ perpendicular to is 
c F X B 

J = -p » 


( 21 ) 


wlierc 

47T F = (div £) E - J E^Ve - j H^'Vy (25) 

Taking the stress tensor of equation 22. as applicable to a collisionless 
plasma (see Appendix) in wliicli 

i= [0, 0, B^(y)], 

E = [0, E (y), EJ, 

— y u 

e = e(y), 
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it is readily shown, as expected, that the electric force is given by 


4-.r ly'-'i caiv n) ly - 4 ry= , 
“ C'-Uv D) 


(29) 

(30) 


Stratton (1941) points out that it is difficult to identify the mechan- 
ical and electric parts of the RUS of equation 21. However in the case of 
plasma tliis appears straightforward as follows. Firstly let us expand the 
right hand side as; 


I 

V 






1 1 1 3b 


-4^ -i CE,. B)] dV, 


(31) 


Then PE - ^ Ej. is the electric force; £ x ^ the Lorentz force due to the 
1 1 . P 

current J; - irVvi is - u V(-jjv), part of the force due to see equation 

gp 

34; x B is the Lorentz force due to the polarisation current, which is 

responsible for accelerating the plasma; ^ x being - £ x £ x £ is an 

additional force exerted by tlie field E on the medium - this force has the 

dimensions, interestingly enough, of (-V • £) £ but is different; notice that 

1 8 

a fox'ce (-£ * £) £ does not appear in this approach, finally (E x B) 

is the rate of increase of the momentum of the electromagnetic field (c.f, 
Chandrasekhar, 1960) . (see discussion later in this paper) « 

It may be noted that equation 21 was derived (Stratton, 1941) under the 
condition that p is a function of position but not of field intensity. This 
is so in all applications in this paper. The value of p in a 2-D collision- 
less plasma (see equation 68a) independent of B for the reason that Pj^/B^ is 
constant as B is changed in such a plasma (Chandrasekhar, 1960) . 

Let us now investigate the steady state flow of plasma as prescribed by 
equations 26, 27 and 28. 
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Case I; Suppose « E^, 

The equations 9 and 24 yield 


By [^(sEy)] - iVW 


3H, 1 


= 0 , 


(32) 


or 


57 CeEy^) - ^ (yH/) - 0. 


(33) 


So 


eE ^ ^ = constant. 

✓ H 


(34) 


Equations 67 and 68 show that in the present special case 


P = 


H 


1 + 


8ttp_l • 


(35) 


So equation 34 yields 


Stt 


8tT 


- p_L = constant 


(56) 


Case II: No restriction on except V x ^ = 0. 

Now equation 24 yields 

''y t57 V - i §7 - T 57 = «■ 

3E 

1 57 ? 5y ^ " “■ t”’ 

This predicts the natural existence of electrostatic fields parallel to 
the magnetic field in plasmas. 
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It follows that 




J. 

3y 


(B^ + 8ttp 


- 


(38) 


To calculate p must first be reduced by V • £, to find div E (equation 11). 
To a high degree of approximation div E. ~ 0, otherwise itolerably large 
electric fields would exist. Therefore, 


BE BE 

^ I 

8z By 


(39) 


Perhaps this is the simplest demonstration of the existence of significant 

8E^ ~ 9E 

electric fields parallel to B in structures in which “-s— ~ 0, and — 5 ^ is non- 
3E 

zero and » The relationships 36 to 39 are testable, in principle, 

with measurements of plasma parameters in essentially laminar plasma flows in 
space above auroras, in the low latitude boundary layer of the magnetosphere 
and in flows in interplanetary space. 

Still with conditions 26, 27 and 28 let us discuss the effect of F 

Zr 

(equation 30) . It is clear that because of the large values of e that occur 
in plasmas (e.g., up to 10 ^ in the magnetosphere and i"iterplanetary space). 
The equation of motion of plasma along B may be v^ritten 


r 

^m dt " 


div D 


— ) E 
4tt z 




Vp , 


(40) 


where g is the component of gravity along the field B_. p is the pressure 
of the plasma measured parallel to £. Clearly the effect of the term fdiv 
£) E^/ 47 T can be significant in structures less than a certain cross section 
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(across B), Usually ;,hen approximate neutrality is assumed in many plasma 

r ft 

systems this term is neglected. However the term is of value EyCg” (s:Ey) + 

r\ 

flows with large variations of in the y direction, 

and correspondingly large contributions are made to F . Examples of this in 

z 

geophysically significant situations are discussed later in this paper. This 
body force would cause jetting of the plasma along B. 


The Current Orthogonal to B 


Expressed generally by equations 25 and 26 the portion of J driven by 


electric field, hereafter called the dielectric current, is given by 

1 


J 1 

-el 


c[(div D) E - i E_[Ve] x B 


47TB2 


(41) 


In the special case given by equations 26, 27 and 28, this current is in the 


x-direction and of magnitude 
c (eE ^3 

j = 


el 


8ttB2 


(42) 


The total current density orthogonal to ^ is, for the same geometry, 

8 r-r- 2 


X 


C ^ - 8ttpx] 

8iP 


(43) 


The Current Parallel to £ 

In addition to the jetting of plasma along ^mentioned earlier the 
parallel electric field will cause current parallel to B_. Such current would 
perturb the simple geometry prescribed by equations 26, 27 and 28. The 
analysis of complete current systems in general magnetic fields is a matter 
for later development. 
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A Numerical Illustration 

Earlier the motions of single charged particles in a homogeneous magnetic 
field of induction ^ and orthogonal electric field ^ which has a gradient VE 
parallel to were studied (Goloj 1976). It was inferred that particles 

with velocity Vj^, given here in c.g.s. units, by 

2 


cB “s 


'Ds 


B > 
s 


(44) 


producing a current in the direction i ^ x £ depending on whether VB has the 

same or opposite direction to (Note that the direction is more properly 

stated as + ^ X B rather than + VE x JB as printed in that paper.) 

“ ^S7E 

Here o3 = Be/m c and ^ = w ^ — tt-tt when q is the charge on particles 

S S S S ITl C 

q 7E B^q^ ^ s 


of mass m„. In the case ,, „ 
s in c 


« g , it is readily shown that 


c^m 


Ds 


2q B3 
‘s 


® 7E2 


cE 

B 


which is charge-dependent and mass-dependent. 

In the special case of a plasma of uniform density this produces a current 


n — V. 
s s c Ds 


eVE^ 

SttB 


(45) 


in agreement with the more general equation 42. 

By comparison of the two approaches, in the case where e is constant, 
equation 42 may be interpreted in terms of a drift for each charged particle 
additional to the well known drift Vj, = c(^ x £) of an amount 

c^m ^E_l^ x ^ 

^VE " 2q B4 • (^^6^ 
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It may be noted that 5 /Vjj is approximately the ratio of the gyroradius 
associated witlt a /olocity to tlio scalo size of the spatial variation of 
Ui . 

J 

In order to ostimato the magnetic change duo to the dielectric current, 
consider equation A2 in the form 

B AB n KyA(eUy) . (47) 

So 

B , 

AB = A (GBy.) . 

Tiie quantity c B /B has the dimensions of velocity V,,, So wo can wi’ito 
y u 

V 

AB : Tp-r B, (48) 

whore is the Alfvdn speed = c/e^. 

Alternatively 

AB ^ V.j^a'irn m B“^. (49) 

Li S S 


where a = 2 if AB/B « 1, and a = 4 if AB ~ B. 

Tabic I sliows estimates of upper limits to AB in various regions of 
interest. 


Region 

Vp(cms“^) 

TABLE I 
B (gauss) 

Mi (amu) 

Ni(cm~5) 

AB (gamma) 

Low lat. Ionosphere* 

105** 

0,3 

20 

5x10° 

4 

Auroral Ionosphere* 

4x105 

0.5 

30 

5x105 

SO 

Mid-Magnetosphere 

3x105 

0.05 

1.5 

1Q3 

0.3 

Low lat. Boundary 

1.5x10^ 

4x10"*' 

1.5 

10 

17 

Layer 

4x10^ 

4X10'** 

1.5 

10 

a»B 

Interplanetary Space 

4x10^ 

B 

1.5 

10 

aiB 

Inner Solar Corona 

10 ^ 

B 

1.5 

10 

ajB 

Bate of Corona 


1 

1.5 

10 ® 

a-B 


* Above altitude at which collision frequency V =* w . 

** Including corotation velocity with earth. s 

12 


ORIGINAL PAGE U 
OF POOR QUALITY 

The effects of dielectric currents are greatest in the boundary layer and in 
interplanetary space and the geophysical consequences of those currents will 
be dealt with in a separate paper. The currents may be detectable in the 
ionosphere particularly in polar regions when the highest spatial variations 
of the quantity eE^ occur. 

It is clear from Table I that the currents must bo taken into account in 
discussions of the physics of the low latitude boundary layer, interplanetary 
medium and the solar corona where the effects of dielectric current can be 
large. 

The current cnr^ integrated in the y direction through the model is, 
from equation 42, 

j(cm~^) = A(BEy^)/87rB emu, (50) 

= n.m. VgV2B. CSl) 

In the case of gradient of electric fields above, or in the vicinity of 
auroras, we assume n. = 5x10® cm m. =30 amu (i.e. NO^), B = 0.5 g, Vp = 
4x10® cm/sec. These are extreme values. From equation 45 then j ~ 1.6 x 
10”** amp cm”^. Integrating over the height of the ionosphere where the ions 
are freely gyrating, i.e., say 180-400 km altitude the total current would be 
3.5 X 10® amp. In the first approximation the altitude variation of the cur- 
rent density is like that of tbe ion mass density. A significant effect 
should hold from 150 km altitude to the altitude of transition of heavy ions 
(NO , 0 ) to light ones (He ,!-!)• Sometimes above auroras this transition 
can be very high in which case the total current could be greater. 

In the low latitude boundary layer (Cole, 1974; Eastman and Hones, 1979) 
in which, e.g., B = 5 x 10"**, = 2 x 10"’ cm~\ n^ = 10 cm”® of H’’’ plus 
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3 cm ^ of He^ the current integrated through the layer is ~ 1.83 x 10 ® 
emu = 1.83 x lO”** amp cm”^. If the boundary layer extends five earth radii 
above and below the equatorial plane, then the total dielectric current flow- 
ing in the boundary layer would be 1.2 x 10® amps. For uncomplicated 
anti- solar flow and velocity gradient pointing outwards, the current is to- 
wards the sun on the dusk side and away from the sun on the dawn side. On 
both sides the Lorentz force of the current is inwards to the magnetosphere. 

The currents can be of magnitude comparable to the conventional Chapman- Ferraro 
currents. Fuller discussions of these problems and applications in the solar 
wind and other astrophysical situations will be taken up in later papers. 

Effects Parallel to ^ 

With the model prescribed by equations 26, 27 and 28, the body force on 
the plasma parallel to ^ is 

The divergence of the displacement vector £ due to spatial variations in the 

y direction together with the large values of e in plasmas create large net 

charge density, acting upon which, E can give a large rate of change of 

z 

momentum to the plasma in the z direction. 

In general we can say that the net acceleration of the plasma is given 
by 

a 

Pmit 

It is of interest to find what length sacle (L*) of variation of 
in the y direction would yield the same acceleration per ionic mass in equa- 
tion 53 as e E It is given by 
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L* == 


“ CB2" 


TA . VpC 
1 E 

Bo 


( 54 ) 


L* is the gyroradius of an ion moving 1 to ^ with a velocity of Vg, Table 
II shows I* for various regions of interest geophysically. 



TABLE 

II 



Regions 

Vg (cms“^3 

B (gauss) 

Mji^(amu) 

L* (cm 

Low lat. Ionosphere 

10® 

0.3 

20 

6.6x10 

Polar Ionosphere 

4x10® 

0,5 

20 

2.6x10 

Mid -Magneto sphere 

3x10® 

0.05 

1.5 

3x10 

Outer Magnetosphere 

Ixio'^ 

0.005 

1.5 

3x10 

Low lat. Boundary Layer 

1.5x10^ 

0.0004 

1.5 

5.6x10 

Interplanetary Space 

4x10^ 

0.00005 

1.5 

1.2x10 


Structures with scale sizes less than these, in the presence of electric 
fields orthogonal to B would necessarily exhibit plasma flows along IB, 
according to equation 53, at places of strong gradient of CeE^,) . This would 
occur naturally at the sharp edges of auroral forms, especially fine auroral 
rays. Such flows may also be induced in the presence of appropriate electric 
fields by chemical releases, e.g. Ba, Li, and discharges from charged-particle 
guns. From 53 it may be inferred that, in a steady state, the distribution 
of plasma along ^ would achieve a distribution given by 

C55) 

The Momentum Per Unit Mass Orthogonal to _B (when pj_ = 0) 

It is clear from the examination of particle motions (Cole, 1976) that 
the current (equation 39) is due to ions in i^roportion to their contribution 
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JN 

to tho mass density of the plasma. When ^ (sEy.^) = 0, the drift of the plasma 
orthogonal to caused by an electric field, is given in a collisionless 
plasma, by 


V 


c E X B 
“'~"p 


C56) 


However in the circumstances given by equations 26, 27 and 28, neglecting 
the momentum of electrons, the momentum per unit mass is given by 


1 8(eE 

V = £. fI3 + X... . 

\ B on 8irdy ^ 


CS7) 


8 

This result means that only if ■gy (eE^^-) « 8 it en E^. may it be stated that, 
approximately. 


cE + VxB = 0 (58) 

In this approximation equation 34 becomes, when e » 1. 

B ^ 1 

Pi •"-§?- 2 Pm'^x' = (59) 

which is a Bernoulli- like relationship, but pertaining to a line orthogonal 
to tho flow. 


The Force - P. x ^ x £ 

This force which was uncovered in equation 31 makes sense as the elec- 
trical counterpart of t.ie force - ^ x ^ x H. Both may be described as forces 
originating from dipoles. The fact that there is no magnetic counterpart of 
(4 it)~^ (div D_3 ^ is related to the lack of magnetic monopoles in the theory. 
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Self-Consistency 

The solf-consistcncy of the method may be verified by calculating the 
forces (electric and magnetic) on the plasma in any uniformly moving frame 
of reference including one moving with the plasma velocity (equation 57) . 

Consider a frame moving with uniform velocity with respect to the 
laboratory (L). Then (Jones, 1964) 

D=Dj^ + VxH/c, (60) 

E = + V X B/c. (61) 

Suppose V « c, so Suppose the plasma is cold i.e., y = 1, so that 

B = il. Also suppose e = constant. For laminar flow such as implied by equa- 
tion 57 


div £ = tliv (62) 

i.e., charge density is the same in both frames. The electric force in a 
frame moving with velocity V multiplied by 4tt is 


(div £) IG = (div £j^) + (div Dj^) V x £/c, 


The current in the frame moving with velocity V is 


div £ 

i = 4 - t— ’ V- 


So the magnetic force in the moving frame is 

div ID 

J X B/c = Jj^ X B/c - ("ifj^) V X B/c 


But 




(63) 


(64) 


(65; 


( 66 ) 
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50 

div D 

E - - 2 X ^/c. 

Consequently there is a balance of forces in all inertial frames. The method 
can claim consistency with the Lorentz invariance of Maxwell's equations. 

The electric currents and forces in a plasma revealed by this method 
come about by allowing the plasma to be driven by external sources and treat- 
ing these sources according to Maxwell's unabridged equations. 

The Magnetic Permeability of a 2-D Collisionless Plasma 

Consider an infinite two dimensional collisionless plasma in which £ = 

0, B^Cy)). In a steady state in the absence of electric field, equation 

51 yields, 

CV x H) X B - i h 2 Vp = 0, C67) 

which in the assumed geometry becomes 

-pH I 7 - I Vp = 0, 
or 

^ CpH") = 0. 

Therefore 

yll^ = K = bVp. 

where K is a constant. 

Also H = h/B. 

In the case of collisionless plasma let us assume that 
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K 5= + Sirpj^. 


This assumption can be proved in the case of low 3 plasma as is done 
below, and of course K is a well known constant of 2-D collisionless plasma. 
Then, 


y = 


1 + 


STTp » 


( 68 ) 


B2 


because p = 1, when p =0. 

JU 

This defines the vectors H and M by 
Slip 

H = (1 + -p- ) 1 = p- 1 


and 


2P, 

M = C- Tt) B = 


4rr i 


The magnetisation vector M is twice the "conventional" magnetisation vector 
and presumably incorporates equal components of the diamagnetism due to 
gyration and due to grad ^ drifting. The total current in the plasma is the 
sum of (477)"^ ^ li c(^ X M), i.e., (477)"^ c V x In the present 

geometry this total plasma current is then 

_c ^ _ c 
477 By " “ B 3y’ 

which is the same as that calculated from particle motions in the field ^ 
(Chandrasekhar, 1560) . 

Consistency with equation (a) is observed because, in this case 
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CV X H) X B = 


K 

47tB 


M = _ J:. 

9y "" " 8 tt 




where 


i = K/B^ 


In the case o£ low 3 plasma, H ~ £ and for present geometry, 

4tt 

Vx H = Vx B = 

Therefore, equation (a) yields 

8pj_ 


1 u2 


9y 8tt 


Vp , 


8tt 9y y- 


Therefore, 


9P1 ” ^ ' 


which leads to the result 


1 STipi 

— = 1 + — nT; — . 

p B2 


(60a) 


It may be of interest to note that this value can be used with the 
classical relationship for the phase velocity of an electromagnetic wave in 
a non-dissipative medium, viz.. 



c^/ep. 


(69) 
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Microscopic plasma physics shows, for a magnetosonic wave, (Stix, 1962) 
moving orthogonal to exactly the same result, i.e., 




1 + 47rp„c^/B^ 
m 

1 +' 3'. 


(70) 


Slip 

where 3, = 

-L n ^ 


Discussion 

The dielectric current orthogonal to a magnetic field in a plasma has 
been put on a general basis. An earlier discussion of the movement of 
charged particles in crossed magnetic and spatially varying electric fields 
revealed a special case of it (Cole, 1976) . The current comes about principally 
on account of the large value of the dielectric constant of a plasma, Struc- 
tures exist in plasmas in the earth’s environment in which the current plays 
a significant role, e.g., in the auroral ionosphere and the low latitude boun- 
dary layer. In a later paper the significance of the current in interplanetary 
space will be demonstrated. This new current should be sought experimentally 
and taken into account theoretically in models of plasma structures. 

The dielectric property of plasma also leads to acceleration of plasma 
bodily along the magnetic field at places of fine structure of the quantity 
and to a new distribution of plasma along 

The coexistence of electric field parallel and perpendicular to ^ in 
plasma flows of special geometry has been demonstrated. Some flows in nature 
appear to have similar geometry. In later papers the consequences of the 
field and currents in the low latitude boundary layer and in interplanetary 
space and the solar corona will be discussed. 
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Those new results are not evident in the widespread approach to plasma 
physics which focusses on internal processes and their description in terms 
of ]E and IB, rendering, for that purpose, £ and H apparently unnecessary as 
vehicles in analyses. However by such examination of internal processes, 
values of e and y can be established. Having these values at our disposal 
allows us to regard a plasma as a medium, with specifiable e and y, which can 
be acted upon from the outside by electromagnetic fields. To analyse such 
action requires the unabridged Maxwell equations employing £, £, £ and H. 

Not, e\?erything which happens in a plasma is the result of an internal 
instability. E.g,, the magnetosphere is an intermediary medium between the 
solar wind "driver" and the ionosphere, sometimes also a "driver". Energy is 
transported across the magnetosphere on time scales far in excess of the time 
it takes for Alfven waves to traverse the system. It is clear that to study 
quasi-steady state macroscopic phenomena such as large scale current systems 
and plasma flows in the magnetosphere, the approach as outlined in this paper 
will need to be taken. The same argument can be made of numerous other systems 
in interplanetary space and stellar atmospheres. 

In the C£j £, £, £) approach to plasmas, all tlie internal currents pro- 
duced by properties of the medium were lumped together in determining y. 

Thus, in the model defined by equations 26, 27 and 28, a value of y was used 
which combined into the vector £ x M both the conventional "magnetisation" 
current of plasma physics and the gradient £ current in a plasma (Chandrasekhar, 
1960) i.e., all the currents related to ion and electron pressure in the 
specified geometry. The value of e is taken from standard plasma physics. 
Parameterising a plasma in this way may facilitate, in some circumstances, 
the macroscopic description of phenomena in it driven from the outside. 
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The virtue of the approach to plasmas in this paper is that it allows an 
internally consistent discussion of non-linear effects in plasmas. However, 
it has relied upon linearised equations (e.g., Stix, 1962) to provide values 
of e and y, There is a sense in which one can speak of a collisionless plasma 
as a perfect Maxwell medium in which the charged particles and the fields, D, 

E, H, mutually interacting determine the behavior of the medium. There 
are no internal strains to be considered as there are in the cases, e.g., of 
solid dielectrics and ferromagnets (see Stratton, 1941) . It is this sense 
which has allowed the interpretation of equation 21 which has been developed 
in this paper. 

It is now possible to see clearly the relationship of the approach of 
this paper to the conventional approach defined by equations 8, 12, 13, 14 
(a and b) and 15. By assuming £ = 0 the conventional approach immediately 
allows only the existence of polarisation current and magnetisation current 
(see equation 12) , The new current defined by equations 24 and 25 is neither 
of these, but is derived from more general considerations. It should be 
appreciated then, that, that approach, though powerful, has its limitations. 

The approach of this paper extends our ability to describe plasma phenomena. 

Generalisation 

Having established the utility of the (^, H) approach, and predicted 

new results with it in special geometry, and observed its consistency with 
other approaches to plasma physics, it is clear that new general results can 
be established by returning to equation 21 and rewriting it as 
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where ijj , is the electric stress tensor and the magnetic one, given by 
ei mag 


4U . = E.D. - i (B 

.• D) 6.., 

(72a) 

^agij = "i“j ■ J ' 

H • B) 6. .. 
~ " 11 

(72b) 

6.. is the Kronecker delta. 






In the most general case assume, using direction cosine notation, that 
E = ECe^^, By, C73) 


and 


8--= BCb^, by. bp 


(74) 


whe 


e^+e^+e^=b^+b^+b^ = l. 
X y z X y c 


Since 


£ = + E| 


C75) 


where 11,1 mean components parallel and perpend icjUar, respectively, to ^ 

D = E[ee,^Cl-b^) + b^e^, ee^Cl-b^) + b^Cy, ee^Cl-b^) + b^e^] (76) 


Then , 


'(’elM = 8? te(l-bp + bp 

- epteCl-bp + bp I 

*elxy = 4f V 


ey^[eCl-bp + by] 


(77) 

(78) 
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while 




i|) d JL e e [e(l-b^J >■ b.J 
^elyx ‘Itt y x x-' 


iJi t3 ^ s' ^rtfl-b ^ + b 1 

Vrtlvvf “ gTT pj- ^yj ^ ^yj 


- c/ [ Cl-bJ + bj 


e/[e(l-b^3 + b^] 


'1^ 1 “ e,o,[e(l>b,) + b„] 

^clvz ‘hr Yz s s 

r * 

'^cizx "" 4¥ ^’x^ 

‘^el.y == 4? ‘^.‘^yteCl-iy . b^J 


el; 


b¥ " V 


e^'[eCl“b,.3 + \] 


I 


- 0^^/ [t'Cl-by) + b^J I 


( 79 ) 
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8ini 


The three relations that are obeyed by a plasma in a steady state are then 


V • W' 1 + ) “ 0 . 

,L '•'♦ol Miuiyy 


U>S) 


Hqiuition 34 is Cound when V x I; h„ ■ b o ^ e^ - 0. Satellite borne 

equipment now measures all tlu'ee components of B and h as well as plasma 
pi'operties from which t and p may be evaluated. So equation S5, using equa** 
tions 77-94 may be tested in steady state structures. The full testing of 
the consistency of measurements of fields B, and and plasma parameters c and 
U aeainst theory for quasi-statie structures would need to allow for the pos- 
sibility of slow but measureablo variations in them. However there would be 
a difficult experimental problem of distinguishing between temporal and 
spatial variations. For this reason it is necessary to have measurements of 
fields and jilasma properties at many different places simultaneously to pro- 
perly understand the pliysics of auroral and ‘’pace plasmas. 

It is of interest to expand equation 71, thus 

fr div D 


nr ^Snech 


0 E t J X i - B_l 


V 

J. 

Sir 


4ir 


B.^Ve 


^ ^ CD X Bij^i V , 


C%) 


and 


. 9E 3M 

^ " it IT IT 

+ 4tt X B - 4itP X V x Fj 

^ 


- -1 

“ at 


ap 


(E X _B ) + 4ir ^ X J1 - 4 it£ x ^ x E_ 


(97) 

(98) 
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The last two terms o£ equations 97 and 98 have bN'jen commented on earlier in 

1 9 

this paper. Chandrasekhar (1962) calls ^ ^ rate o£ change o£ 

1 9 

electromagnetic momentum. This term is a combination o£ 4 ^ CE_ x H), 
which is the Abraham and von Laue expression £or rate o£ change o£ electro- 
magnetic momentum in material media (Stratton, 1941) , plus £orces on the 

medium due to the interaction o£ the electric £ield (^) and the magnetisation 

1 8 

£ield (M) o£ amount ~ ^ ^ corresponds to the interpretation o£ 

£ X H , as the Poynting vector. In slow time variations as envisaged in 

1 9 

this paper, this author considers that ~ ^ (^ x M) should be interpreted as 

a £orce on the medium contributing to its acceleration. 

For variations in the plasma, slow compared to the ion gyroperiod, the 
19 19 

£orces — ^ (E x M) and — ■ 5 ^ (P x B) may be considered to be equivalent in 
C aX. — — C dt “ “ 


their mechanical e££ects to a current in the plasma o£ amount given by 
i* = Cl ^ M ^ X B)] X B B“^. 


(99) 


This expression contains, as one component, the well-known polarisation cur- 
rent. The three new currents implied by equation 99 are 

9E 

j_i* = B"2 B X (M X ^), (100) 

8 M 

j_ 2 * = B“"(E X -g^) X B, (101) 

j_ 3 * = c B"^ B X (P X ( V X E)) . (102) 

9D 

The whole o£ j_* would be in -g:£- and J o£ equation 2. The absence o£ the 
current o£ equations 100 , 101 and 102 in much conventional plasma physics 
which arbitrarily puts £. = 0 should give us pause to consider the assumptions 
o£ that method. Of course ji* and j 2 * are zero in a cold plasma while j_ 3 * 
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is ^ times the well-known polarisation current where A denotes pertur- 

bation (by a wave, say) . 

Now 


cP ^ ~ ^ M 
B ■ AP ~ B ’ AE 




ph 


where is the phase velocity o£ the wave. So, in many circumstances ^ 3 * 
will be a negligible current, but not always. Likewise, if we compare ji* to 
c V X M the ratio is Vg AE/(EC^), again often negligible, but not always. 

As regards the ratio of jz* to c V x M, it is Vg AM/ (MC^) which is 

often negligible. 

Accepting the Abraham and von Laue expression for electromagnetic 
momentum, even in a material medium (Stratton, 1941), the rate of change of 
mechanical momentum per unit volume of plasma would be given by the equation. 


^ (p V) = pE + - J X B - E^ ,Ve - ^ H=’- Vy 

dt — C — — 8 tT JL — oTT — 

1 BP ^ 1 1 

+ — 7rrxB+— PX7TT- + — Ex-^-r-— Mx^ 

c 9t — c — 3t c — 9t c — 9t 


(103) 


The conclusions expressed in equations 99-102, based upon the Abraham 
and von Laue expression for electromagnetic momentum would be important in 
various circumstances, e.g., the current j_ 3 * could be significant in col- 
lisionless shocks, and the currents ji* and j 2 * in relativistic plasmas. 


It should be noted that equation 95 is a "weak" condition for a steady 
state because it does allow the exchange of momentum between mechanical and 
electromagnetic components within the volume considered. From equation 103 
and equation 2 it can be inferred that a "mechanical" steady state of the 
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plasma is acquired when the following stronger conditions exists 


(div D ) E + C V X H) X B - i Ej^^Ve - j H^Vp = 0 


(104) 


For the conditions prescribed by equations 73, 74 and 75, equation 104 yields 




1 , a 1 3(B /y) 3CB /|i) 

* T al \ ^ > 


3z 


+ B . 

y 9y 


9(B /y) 9(B/y) 

^ 1 


9x 


C105) 


and similar equations for the y and z components of 104. It is clear that 
these can be integrated only for the simplest of geometry such as that of 
equations 26, 27 and 28. From an experimentalist's viewpoint this would mean 
that only with the simplest of geometry may the basic equations be tested in 
the steady state, in integral form, using single point measurements. The 
alternative is to test them in differential form (like equations 104) . In 
space physics terms this calls for clusters of satellites or space probes 
identically instrumented to study space plasmas definitively. 

There is much discussion in the literature as to the appropriateness of 

E X H £ ^ i . , 

■ (the Abraham and von Laue expression ) or - (the Minkowski expres- 

sion) for the electromagnetic momentum in a material medium (see Mjziller, 
1972). Miller suggests that it may not be possible to resolve this question. 
The question may be more than one of the semantics of equivalent descriptions 
of forces on a system consisting of electromagnetic fields and matter. In 
the case of plasma in a magnetic field the so-called "polarisation current" 

is clearly seen to play a role in accelerating plasma (c.f. , Spitzer, 1962). 

1 

But the force due to the "polarisation current" is ~ x ^ and this is part 
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of -r— •?— CD X B1 , This indicates that when we come to discuss the mechanical 
momentum of a plasma the semantic dilemma is broken and that the Abraham and 
von Lauc expression for electromagnetic momentum is valid. This points to the 
reality of the new currents implied by equation 99. Moreover this discussion 
points to the value of space plasma measurements in answering fundamental 
questions in electromagnetism, because in collisionless plasmas, not only can 
the fields and particle distributions now be measured, but this medium has 
the virtue of the lack of additional internal forces to contend with. What 
is lacking at present is measurements of such plasma with sufficient spatial 
and temporal resolution. 

Finally it should be mentioned that it is expected that the new current 
(equation 42) should also exist in collision dominated plasma because a 
gradient of electric field energy density is still equivalent to a force on 
the plasma. Further, it should be understood that the currents uncovered in 
this paper must imply that stability theory of plasma is far from complete 
and these currents should be taken into account in such theory. 
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The stress tensor 

As Stratton 0941) explains, the stress tensors would be different from 
those of equations 72a and 79b, if e were a function of E, and p a function of 
B. since 


4ir n m c 
. . ^ s s 
E = 1 + ^ , 

n 


e is not a function of E and the electric field energy density is given by 


u= / 


E » d p 
8tt 


E 


= e 


EdE 

Btt 


e E 
8tt 


Also since 


P 


StrPi 

1 + — — = 1 + p , y is not a function of B for 2- 


mensional adiabatic changes, because, in such changes, n/B is constant 
and Tj^/b is constant ( chandrasehhar 1960). 

It follows that the magnetic energy density is given by 


w 


B 

/ H 

o 


, B 2 

- f B • dB = . 

y ^ - - y 


Therefore, it is appropriate to employ the stress tensors 72a and 72b in the 
case of 2-D collisionless plasmas. 
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